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SUMMARY 

The onset and groirth of delaminations in unnotched [±30/±30/90/90] 

s 

graphite-epoxy laminates is described quantitatively. These laminates, 

designed to delaminate at the edges under tensile loads, were tested and 

analyzed. Delamination growth and stiffness loss were monitored nondestruc- 

tively. Laminate stiffness decreased linearly with delamination size. The 

strain energy release rate, G, associated with delamination growth, was 

calculated from two analyses. A critical G for delamination onset was 

determined, and then was used to predict the onset of delaminations in 

[+45 /-4-5 /O /90 ] (n=l,2,3) laminates. A delamination resistance curve 

n n n ns 

(R-curve) was developed to characterize the observed stable delamination 
growth under quasi-static loading. A power law correlation between G and 
delamination growth rates in fatigue was established. 





INTRODUCTION 


A commonly observed failure mode in laminated composite materials is 
delamination between the composite layers. Delaminations may develop during 
manufacture due to incomplete curing or the introduction of a foreign 
particle; they may result from impact damage; or they may result from the 
interlaminar stresses that develop at stress-free edges or discontinuities. 
Furthermore, delaminations may grow under cyclic loading. Delamination 
growth redistributes the stresses in the plies of a laminate, and may 
influence residual stiffness, residual strength, and fatigue life. Hence, 
a fatigue analysis for composite materials should take into account the 
presence and growth of delarainations . 

A bibliography of experimental and analytical work on delamination is 
contained in reference [l]. One of the most promising techniques for 
characterizing delamination growth is based on the rate of strain energy 
released, G, with delamination growth [l]. Previous work [2] has shown 
that the cyclic growth rate of debonds between the metal and composite 
components of reinforced panels could be correlated with G. Measured 
critical G values have been used in sophisticated analyses [3,^] to pre- 
dict the onset of edge delaminations in unnotched composite laminates. 

In the present study, a simple technique was developed, employing 
strain energy release rates to characterize the onset and growth of 
delaminations in a composite laminate. First, the damage that developed 
in unnotched [±30/±30/90/90] graphite-epoxy laminates under static 
tension loading and tension-tension fatigue loading was determined. Next, 
stress distributions generated from a finite element analysis were 
correlated with the observed damage. Then, during quasi-static test 



loads, delaraination growth and stiffness loss were monitored nondestruc- 
tively to relate laminate stiffness and delamination size. The resulting 
test data and analysis were used to derive a closed-form equation for the 
strain energy release rate, G, associated with delamination growth. Next, 
a critical value of G for delamination onset was determined. It was then 
used to predict the onset of delamination in 

laminates. A delamination resistance curve (R-curve) was developed to 
characterize the observed stable delamination growth during quasir-static 
loading. Finally, a power law correlation between G and delamination 
growth rates in fatigue was established. 

SYMBOLS 

[a] extensional stiffness matrix 

A delaminated area 




V 


A* 

a 

Aa 

da 

dN 

[B] 

b 

C,3 

[D] 


area of interface containing a delamination 
strip delamination size 
incremental strip delamination size 
delamination growth rate in fatigue 
coupling stiffness matrix 
half-width of laminate cross section 


empirically determined coefficients 
bending stiffness matrix 
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E 


axial stiffness of a partially delaminated laminate 


®LAm axial laminate stiffness calculated from laminated plate theory 

E* axial stiffness of a laminate completely delaminated along one or 

more interfaces 

Ej’ axial stiffness of the i"*^^ sublaminate formed by a delamination 

Eq initial tangent modulus of an undamaged laminate 

%1»^22’^33 moduli 

^ rate of stiffness change with delamination area 

dA 

— rate of stiffness change with fatigue cycles 

dN 

^11 5*^12 ’^13 shear moduli 

G strain energy release rate associated with delamination growth 

s'train energy release rate components due to opening, 
in-plane shear, and out of plane shear fracture modes 

GjpgM values of Gj, calculated from finite element 

analysis 

G^ critical strain energy release rate for delamination onset 

Gjj delamination resistance 

maximum strain energy release rate in constant strain 

iiiax 

amplitude fatigue test 
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U,V,¥ 
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du 

dA 

M 

dA 

x,y,z 

X,Y,Z 

^11 


e 


ply thickness 

saturation spacing of cracks in 90° plies 
gage length used to measure axial displacements 
number of fatigue cycles 
number of plies in a laminate 
laminate thickness 

thickness of i^^ sublaminate formed by delamination 
displacements in x,y,z directions 
displacement functions in x,y,z directions 
material volume 

rate of strain energy released as flaw extends 
rate of work done by applied load as flaw extends 
cartesian coordinates 
nodal forces in x,y,z directions 

first element of the inverse extensional stiffness matrix 
(i,J = 1 , 2 , 3 ) 

nominal axial strain 

uniform axial strain assumed in finite element analysis 
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nominal axial strain at onset of delamination 
^max maximum cyclic strain level in fatigue r 

^12’^13»^23 Poisson's ratio * 

c?x axial stress in a ply 

remote axial stress applied at onset of delamination 
interlaminar normal stress between plies 
SPECIMENS AND APPARATUS 

Unnotched [i30/±30/90/90]^ , T300-5208 graphite-epoxy laminates were 
tested in tension. This laminate was designed to have relatively high 
tensile interlaminar normal stresses at the edges resulting in the formation 
of a delamination [5], Figure 1(a) shows the delamination that developed 
along the edge. 

The specimens were 254 mm (lO in.) long by 38 mm (1.5 in.) wide. These 
eleven-ply laminates had an average ply thickness of 0.l4 mm (.0054 in.). 

Specimens were tested in a closed-loop hydraulic testing machine. The 
specimen length between the grips was I80 mm (j.O in.). 

Also shown in figure l(a) is a pair of linear variable differential 
transducers (LVDT's) , that were mounted on the specimen to measure displace- 
ments over a 102 mm (4.0 in.) gage length. To prevent slippage, a fast 
drying glue was applied to the central 5.6 mm (.22 in.) portion of the 

LVDT mounts where they touched the specimen. Dviring "strain-controlled" f 

loading, these LVDT's were used as the feedback device in the closed loop. 

r 
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Dye-penetrant enlianced radiography was used to monitor delamination 
growth through the specimen width. Diiodohutane (DIB), a dye penetrant 
opaque to x-rays, was injected along the delaminated edge. The film was 
placed immediately behind the specimen. While still mounted in the test 
machine, specimens were exposed to x-rays generated for five seconds at l8 kV 
from a portable point-source \init positioned 386 mm (15«5 in.) away from the 
specimen. The radiographs showed the location of the delamination front 
(fig. 1(b)). The dark outline in the center of the picture is the "shadow" 
(x-ray image) of the LVDT rods used to measure displacements. 

DAMAGE DEVELOPMENT 

The same type of damage developed during both quasi-static tension 
and constant-amplitude, tension-tension fatigue. First, a few isolated 
cracks formed in the 90° plies. These were followed almost immediately by 
the formation of small delaminations along the edge, as seen in figure 2(a). 
When the delaminations formed, the number of 90°-ply cracks increased 
significantly along the delaminated length of the specimen. Most of the 
90°-ply cracks , which appear as horizontal lines on the radiographs in fig- 
ure 2, extended beyond the delamination front outlined in the x-ray photo- 
graph, Many of the ply-cracks immediately extended halfway across the 
specimen width. As loading continued, additional delaminations formed and 
joined with original delaminations. Delaminations grew much more rapidly 
along the length of the specimen than across the width (fig. 2(b)). 
Eventually, two delaminations, one on each side of the specimen, extended 
along the entire specimen length between the grips (fig. 2(c)), after which 
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tlie deXaBiinations continued to grow across the width. Loading was terminated 
when the delamination front reached the shadow of the LVDT rods (fig. 2(d)). 

To illustrate the location of damage through the thickness, a few 
acetate tape replicas of a delaminated edge were made [6], Figure 3 shows 
photogra,phs of two replicas and a portion of the delaminated edge. As shown 
in figure 3(a), ply cracks extended through the thickness of all three 
interior 90 plies. As shown in figure 3(h), delaminations formed and grew 
in -30O/90° interfaces, typically shifting from one interface, through 
90° ply cracks, to its symmetric -30°/90° counterpart. However, delamina- 
tions did not shift interfaces at every 90° ply crack encountered. As shown 
in a photograph of the delaminated edge, figure 3(c), interface shifting did 
not occur in a regular pattern. 

Besides formation of delaminations at the -30/90 interfaces and cracks 
in the 90° plies, an occasional angular crack formed in the innermost -30° ply. 
The replica of the edge in figure 4(a) shows two angular cracks in the inner- 
most -30° plies originating at a 90° ply crack tip and creating delaminations 
in the +30/-30 interfaces. As shown in the radiograph in figure 4(b), the 
delaminations in the +30/— 30 interfaces were small and triangular in shape. 
These +30/-30 delaminations often temporarily arrested initial -30/90 inter- 
face delamination growth along the length of the edge. However, -30/90 inter- 
face delaminations eventually Joined up and grew, whereas isolated 
+30/-30 delaminations usually remained small. 

STRESS ANALYSIS 

Two approximate analyses were, used to obtain quantitative predictions 
of the onset and growth of delaminations. The first was a quasi-three- 
dimensional stress analysis that yielded stress distributions and strain 



energy release rates. The second was a simple rule of mixtures analysis 
that was used along with laminated plate theory to calculate stiffness loss 
and strain energy release rates. 

Because delaminations form in xinnotched laminates as a result of the 
interlaminar stresses that develop at the edge, a q.uasi-three-dimensional 
finite element analysis [7] was performed. The finite element analysis was 
used to calculate stress distributions in the [±30/±30/90/^]_ laminate 

o 

for a unit axial nominal strain (e^ = l). Some details of the analysis are 
described in Appendix A. 

Figure 5 shows that the through-thickness distribution of the inter- 
laminar normal stress, calcvilated at the edge is compressive in the 

outer 30° plies but reaches a relatively high tensile value at the -30/90 
interface and throughout the 90° plies. Also shown in figure 5 is the 
approximate distribution through the thickness calculated from 

laminated plate theory and an assvtmed stress distribution across the 
width [8]. This plot also shows the highest tensile 0 „ stresses to be at 
the -30/90 interface and within the 90° plies. 

Figure 6 shows a distribution of a across the specimen width, near 
the edge, at the -30/90 interface, as well as a distribution of the axial 
stress, 0^, in the adjacent 90° ply. Both and have high tensile 

values at the edge. 

These stress distributions showed reasonable correlation with the 
observed damage that developed. Indeed, examining and interlaminar 

shear stress distributions are helpful in identifying likely delamination 
sites. However, interlaminar stress distributions calculated from finite 



element analyses were not usefxxl for modeling damage growth quantitatively 
because the magnitude of calculated peak stresses at the edge varied with 
mesh size. Furthermore, linear elastic analysis suggests that the inter- 
laminar stresses at ply interfaces can become singular at the edge [T]» This 
singular behavior would preclude the use of a failure criterion based on 
maximum interlaminar stress values. Therefore, an alternate approach, based 
upon strain energy release rates, was adopted to quantitatively describe the 
onset and growth of delaminations. 

STIFFNESS LOSS 

In many composite laminates, stiffness loss may reflect delamination 
growth. Furthermore, the rate of stiffness loss with delamination growth can 
be directly related to strain energy release rates. Therefore, analysis and 
experiments were performed to correlate laminate stiffness and delamination 

size. 

Rule of Mixtures Analysis 

To analyze stiffness loss due to delamination, a simple rule of mixtures 
analyses, along with laminated plate theory, was used. First, the stiffness 
(tangent modulus) of a balanced, symmetric composite laminate (fig. Ta) was 
calculated from laminate theory [9 >10] as 


where first element of the inverse extensional stiffness 

matrix, (i,j = 1,2,3), and t is the laminate thickness. Next, 
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assuming a complete delamination in one or more interfaces, and using the 
rule of mixtures assumption that the suhlaminates formed' undergo the same 
axial strain (hut no longer have the same transverse strains), results in 


E* = 


n 


i=l 


t 


( 2 ) 


where 


E* - stiffness of a laminate completely delaminated along one or 
more interfaces 


Ej_ = the laminate stiffness of the i"^^ sublaminate formed by the 
delamination 


t^ = the thickness of the i't^ sublaminate 

Although eciuation (2) represents a two-dimensional formulation, E* will 
depend upon which interfaces delaminate. This, in turn, determines the 
stiffness, E^, and thickness, tj^, of each new sublaminate. Hence, eq.ua- 
tion (2) is sensitive to the through-thickness location of the delamination. 
For the [±30/±30/90/90]g laminate, assuming a delamination in both -30/90 
interfaces (fig. 7b) equation (2) becomes 

E. . 3E(9o)^ _ 

• 11 

% Finally an equation for the stiffness, E, of a partially delamina,ted 

specimen was developed using the rule of mixtures. Equal— sized delaminated 
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strips were assumed to exist at both edges of the lamiinate (fig. Tc). Then, 
by assuming the laminated and delaminated portions of the specimen act as 
independent components loaded in parallel, the rule of mixtures yields 

® - ®lam) ^ w 

A more general form of equation (4) may be developed by assuming that the 
relationship between laminate stiffness loss and delamination size can be 
represented by 

^ ~ ^LAM = ^ /,-N 

E« - Elam a* 

where 

A = delaminated area 

A* = total interfacial area 

Rearranging equation (5) yields 

E = (e* - Ej^^) + Ej^ (6) 

Equation (4) is a special case of equation (6) where 
S. = iL- 

b A* 


f 


f 
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Experiments 


To verify the linear relationship between stiffness and delamination 
size implied by equations (4) and (6), four quasi-static tensile tests were 
conducted. The specimens were loaded in a strain controlled mode until a 
delamination formed. Then, the specimens were unloaded to ten percent of 
the peak nominal strain, DIB was placed on the specimen edges, and an x-ray 
photograph was taken. Next, the specimens were reloaded in increments 
of 250 ym/m above the previous maximum strain level. This procedure was 
repeated until the specimen was almost totally delaminated. 

During each loading, output signals of the two LVDT's were averaged, 
and load deflection curves were plotted on an X-Y plotter. The initial 
linear portion of each plot was used to calculate laminate stiffness 
corresponding to the damage recorded in the previous x-ray photograph 
(fig. 8). Delaminated areas recorded on the photographs within the 102 mm 
(4 in.) gage length were measured with a planimeter. To minimize data reduc 
tion error, each delamination was traced three times and measured areas were 
averaged. Then, a strip delamination size, a, having equal area over the 
LVDT gage length as the measured delamination, was calculated (fig. 9). 

Figure 10 shows a plot of normalized stiffness, E/E^, as a function of 

normalized delamination size a/b. A least-squares regression line for the 

data indicated that E* = 0 .t42Eq, where E^ is the initial tangent modulus 

measured. Hence, a total delamination in the [±30/±30/90/90]„ laminate 

s 

would result in a 25.8 percent reduction in laminate stiffness. The data 
agreed with the linear rule of mixtures equation (4), normalized by Ej^^, 
where E* was calculated from equation (3) and sublaminate stiffnesses were 
calculated using equation (l) (see appendix B) with material properties 
from reference [ll]. 
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Finally, eq.uations (2) and (H) were used iteratively to calculate 

laminate stiffness for specimens having -30/90 interface delaminations and 

concurrent, although small, +30/-30 interface delaminations (see Appendix C) 

The contribution of 90° ply cracks to laminate stiffness loss was also 

considered (see Appendix C). However, the net effect of both secondary 

mechanisms (+30/-30 delaminations and ply cracks) on stiffness loss was 

negligible for the [±30/±30/90/90]„ laminate. 

s 

STRAIN ENERGY RELEASE RATE 

For an elastic body containing a planar flaw of area A, the strain 
energy release rate, G, is the difference between the rate of work done, 
dW/dA, and the rate at which elastic strain energy is stored, dU./dA, as the 
flaw area increases [l2], i.e.. 


Q _ dW _ du 
dA dA 


( 7 ) 


Assuming that a nominal strain, e, is sufficient to extend the flaw, the 
work term vanishes. Then, if u is expressed as a product of the strain- 
energy-density and volume of the body, (/, substituting Hooke's law into 
equation (j) yields 


G = -1/ ^ 

2 dA 


( 8 ) 


where 


dA 


is the rate of stiffness change as the flaw extends. 



In this study, the body was a tensile-loaded, unnotched composite 
laminate, containing edge delaminations. The strain energy release rate 
associated with the growth of edge delaminations can be calculated by 
assuming two strip delaminations (fig. 9) where 

1/ = 2b il t 
A* = 2b A 
k = 2 % a. 

dA = 2£ da 


Then, substituting equations (9) into equation (8) and differentiating 
equation ( 4 ) yields 


G = /E 

2 LAM 



( 10 ) 


Equation (lO) may also be derived for an arbitrary-shaped delamination by 
substituting equation (6) into equation (8), differentiating, and noting 
that 1/ = A*t , Hence, as indicated in equation (lO), the strain energy 
release rate associated with delamination growth is independent of the 
delamination size. The magnitude of G depends only on the laminate layup 
and location of the delaminated interface(s) (which determine and E*), 

the nominal strain, e, and the laminate thickness, t. 

Furthermore, the strain energy release rate (eq. (lO)) may have con- 
tributions from any of the three components ^III’ corresponding 

to the opening, in— plane shear, and out-of— plane shear fracture modes. In 
addition, near the edge, G may deviate from the value predicted by 
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equation (lO), which was developed using laminated plate theory and the rule 
of mixtures, Therefore, a virtual crack extension technique was used with 
the quasi-three-dimensional finite element analysis to calculate Gj, 

and Gjii as a function of delamination size (see Appendix A). As shown / 

in figure 11, the finite element analysis indicated that the total G, 
represented hy Gj plus Gjj negligible), reached the value pre- 

dicted from equation (lO) once the delamination had grown a very small dis- 
tance in from the edge. 

Because G increased rapidly with "a" near the edge, a small delamina- 
tion that formed (for whatever reason) at the edge would be expected to 
undergo rapid initial growth. This behavior was observed in the quasi-static 
tension tests used to generate stiffness data. As soon as a delaraination 
was detected, the loading was stopped. The formation and growth of the 
delamination to some finite size appeared to be nearly instantaneous. There- 
fore, G calcilLated from equation (lO) at the nominal strain where delamina- 
tion was first detected was considered to be the critical value, G^ , required 
to form the delamination. This G^ was then used to predict the onset of 
delamination in other laminates. 

\ 

DELAMINATION ONSET 

To predict the onset of delamination in other laminates, several things 
were done. First, tension tests of [±30/±30/90/90] laminates were run to 
determine the nominal strain level, at which delamination begins. 

Next, was used in equation (lO) to predict a critical G^ for the onset 

ir 

of delaminatiorio Then, was used to predict the nominal strain at the 

onset of delamination in other laminates. A more detailed description of / 

the procedure follows. 
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Critical Gq Determination 


First, eighteen [±30/±30/90/90] graphite epoxy laminates were loaded 

o 

monotonically in tension at a rate of 4^1.5 N/sec (lO Ibs/sec) until a 

delamination was detected. The load level corresponding to delamination onset 

was recorded and the corresponding applied stress, a^, was calculated. Then, 

to determine the nominal strain at the onset of delamination, e , a_ was 

c c 

divided hy calculated from laminated plate theory using the following 

elastic properties from reference [ll] 

Ell = 138 GPa (20.0 Msi) 

Egg = 15 GPa (2.1 Msi) 

G;Lg = 5.9 GPa (0.85 Msi) 

Vj^g “ 0 . 21 

The average e value was 3^70 ym/m. In addition, e was determined from 
c c 

LVDT measurements on the four tests conducted to generate stiffness data. 

In each of these four tests, the load deflection plot was linear until the 
delamination formed. The average value of where the load-deflection 

curve deviated from linear for these four tests was also SkfO ym/m. 

Next, e was substituted into equation (lO) to determine G„. Stiff- 
ness were calculated from equations (l) and (3), respectively, 

using elastic properties from reference [ll]. The average laminate thick- 
ness, measured with micrometers, was I. 5 I mm (0,059^ in.). A value 
of 137 J/m^, ( 0.78 in-lbs/in^) was calculated for the critical strain energy 
release rate. 
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Delamination Onset Prediction 


To predict the onset of delamination in other laminates, eq.\iation (lO) 
was inverted to yield 

The critical detennined from [±30/±30/90/90] laminate data was used 

c s 

in equation (ll) to predict the nominal strain at the onset of delamination 
in “ 1>2,3) T300-5208 graphite-epoxy laminates 

having the same stacking sequence but different thicknesses. To evaluate E*, 
delaminations were modeled along both 0/90 interfaces where high tensile a 

z 

stresses were anticipated [3] on the basis of the approximate analysis of 
reference [8]. Delaminations have been observed in the 0/90 interface of 
a [±45/0/90]^ laminate in reference [l2]. Then, stiffnesses and E* 

were calculated from equations (l) and (2), respectively, using elastic 
properties from reference [ll]. / 

A ply thickness of 0.15 mm (0.0057 in.) was determined in reference [3] 
for the 8-ply (n = l), l6-ply (n = 2), and 2l|-ply (n = 3) laminates. In 
figure (12), predictions of were compared to values calculated by 

dividing measured values from reference [3] by Ej^. Because the data 

from reference [3] represents the average of only two or three tests for 
each stacking sequence, only preliminary conclusions can be drawn. Never- 
theless, the good comparison indicates that may be independent of the 

ply orientations that make up the delaminating interface. For example, 

Gq may be the same for delamination onset in the -30/90 interfaces of 
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[±30/±30/90/90]„ laminates and for delamination onset in the 0/90 inter- 
faces of the [+^5^/-^5n/On/90^]g (n = 1,2,3) laminates. In fact, equa- 
tion ( 11 ) indicates that the dependence of on stacking sequence and the 

ply orientations that make up the delaminating interface is accounted for 
in the - E*) terra. 

In addition, both equation (ll) and the data indicate that delaminations 
will form at a lower nominal strain in thicker laminates of identical stacking 
sequences. According to reference [4] interlaminar stress distributions 
calculated from elastic analysis will be identical for all [+45„/-45„/0„/90„ ]„ 
laminates subjected to the same applied stress. Therefore, the authors con- 
cluded that a failure criterion, if based on critical interlaminar stresses, 
would not predict the thickness dependence of delamination onset. Hence, 
the strain energy release rate appears to be the most useful parameter for 
quantitatively predicing the onset of delamination after the delamination- 
prone interface(s) have been identified from a stress analysis. 

DELAMINATION GROWTH 

Quasi-Static Tension 

Edge delamination has been observed to be a stable fracture process 
in laminates subjected to tension loading [l,5,l4]. Hence, the applied 
load must be increased to force the delamination to grow. Tension tests on 
the [± 30 / 130 / 90 / 90 ]^ laminates confirmed this observation. Unstable growth 
of the delamination through the width did not occur before the laminate 
failed (fractured into two pieces). In the four quasi-static tests con- 
ducted to generate stiffness data, delaminations did not grow after the mean 
S'PPlisd load exceeded 13,350 N (3000 lbs). The four laminates eventually 



failed at a mean load of 20,000 N (1*500 lbs). Therefore, the 

[ 130/130/90/90 ] g laminate is well suited to studying stable delsimination 
growth . 

The stable growth of flaws can be characterized using the crack growth 
resistance curve (R-curve) concept of fracture mechanics [15], Therefore, 
a delamination resistance curve was constructed. Strain energy release 
rates, G, were calctdated from equation (lO). Because G does not depend 
on delamination size, it appears as a horizontal line in figure 13. The 
three horizontal lines shown are G values calculated for a single specimen 
at three successive nominal strain levels. In addition, the delamination 
resistance, Gj^, was calculated using the maximum nominal strain in equa- 
tion (10). However, Gj^ was plotted as a point corresponding to the size 
(see fig. 9) of the delamination created by the nominal strain. As shown 
in figure 13, the curve formed by all such points, generated during the 
four [±30/±30/90/90]g quasi-static tension tests, constitutes the delamina- 
tion growth resistance ciirve (R— curve) for the graphite epoxy specimens 
tested. The critical G^ used in the previous section to predict the onset 
of delamination represents the first value of G^, i.e., the first point on 
the R-curve. Hence, the R-curve characterizes the laminates resistance to 
delamination growth under tensile loading. 

If the R-curve is independent of the ply orientations that make up 
the delaminated interface, like G^ appears to be (see previous section), 
then equation (lO) and the R-curve can be used to predict the growth of 
delaminations under quasi-static tension in other laminates. As the 
loading is increased, G can be calculated from equation (lO) and 



compared with the R-curve to predict delamination size. Hence, delamination 
size as a function of applied load could be predicted for other laminates. 

Fatigue 

ExpzAlrmnti) 

Constant amplitude, tension-tension, strain-controlled fatigue tests 

of [±30/130/90/90] graphite epoxy specimens were conducted at a frequency 

s 

of 10 Hertz and a strain ratio of 0.2, Specimens were loaded slowly in 
tension until a small delamination appeared. Laminate stiffness was measured 
during this initial loading. Then, the specimen was tinloaded and delamina- 
tion size was recorded using DIB-enhanced x-ray photography. Next, the 
specimen was reloaded to the mean strain and the stiffness of the delaminated 
specimen was recorded. Then, constant-strain-amplitude cyclic loading was 
applied. The cyclic loading was interrupted at specified intervals to 
measure delamination size and static stiffness. Hence, a photographic 
record of delamination size and a record of static stiffness as a function of 
load cycles were acciomulated. Cyclic loading was terminated when the 
delamination had grown across most of the specimen width and the front was 
obscured by the x-ray image of the LVDT rods (fig. 2 (d)). 

V&lamination cjAowtk mtz mea6uA.m&nt6 

Figure l4 shows a typical plot of delamination size as a function of 
load cycles. A strip delamination size, a, was calculated from the delami- 
nated area (fig. 9) measured from x-ray photographs using a planimeter. 

Three separate planimeter tracings of each delaminated area were performed to 
minimize data reduction error. As indicated in figure l4, once the delamina- 
tion had grown over the entire length of the specimen edge, a constant growth 



rate, da/dN, was eventually achieved. Table 1 lists the growth rates, 
da/dN, determined from least squares linear regression analysis of the data 
for each fatigue test. 

Figure 15 shows- a typical plot of static stiffness loss as a function of 
load cycles. This plot also became fairly linear once the delamination had 
grown away from the edge. Stiffness degradation rates, dE/dN, were calcu- 
lated from least squares linear regression analysis of the data over the 
same cyclic range used to fit da/dN data. Because stiffness was found 
to be linearly related to delamination size, differentiating both sides of 
equation (4) with respect to n-umber of cycles (N) allowed an alternate 
determination of delamination growth rates, da/dN, from measured stiffness 
degradation rates, dE/dN, i.e.. 


da _ / b \ dE 

as “ I®* - Slam] a“ ’ 


(12) 


Therefore, delamination growth rates for the edge delamination specimens 

could he estimated without measuring the delamination size directly (Table l). 

Hence, the [±30/±30/90/90 ] specimens used in this investigation should 

s 

be useful for generating baseline delamination growth data* 

Vcuta ao/iAoZcutton u)XXh anali^6tii 

Four fatigue tests were conducted at each of three maximum cyclic 

strain levels. The maximum cyclic strain levels, chosen were 3000, 

3250 , and 3500 ym/m. The maximum strain energy release rate, G , was 

luax 

calc’ulated from equation (lO) using Therefore, a constant 

was also a constant G^^^ test. The mean of the four delamination growth 

max 



rates, da/dN, was determined for each of the three G levels. A power 
curve of the form — = c G^^^ ^ was fit to the three mean values of da/dN 
and the three values of G^^^ using a least squares routine. An excellent 
correlation was achieved for both delamination growth rate measurement tech- 
niques . 

Figure 16 shows — as a function of G^^ ^ for growth rates estimated 
directly from measured delaminated areas . The growth rates for all four 
tests at each level are shown along with the mean growth rate. Fig- 

ure 17 compares the least-squares power law fits using direct area measure- 
ments and indirect stiffness estimates of delamination growth rates. Because 
both da/dN and G can be calculated without direct measurement of delamina- 
tion size, the stiffness technique presents a relatively simple means of 
generating data to determine the dependence of empirical parameters c and 3 
on differences in load history, frequency, temperature, etc. 

With the parameters c and 3 determined, the power law might be 
to other laminates. The calculated strain energy release rates 
could be used with c and 3 measured from edge-delamination baseline 
tests to determine delamination growth rates. 

CONCLUDING REMARKS 

Methodology Summary 

A methodology for analyzing the onset and growth of delaminations in 
composite materials was formulated based on findings in the current work 
and in the literature. The methodology is as follows. 

First, a stress analysis of the particular material, configuration, 
and loading must be performed. The analysis should establish where 
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delaminations will be located. For 'unnotcHed laminates, the approximate analyses 
discussed in this paper [8,l6] might be adequate. After the location is established, 
the delamination onset and growth can be characterized quantitatively Using 
strain energy release rates. Such a characterization incorporates the influ- 
ence of material volume. Furthermore, determination of the singular stress 
field at the delamination front is not required. 

The current work indicates that a total strain energy release rate 
may be sufficient to characterize the onset and growth of edge delaminations 
in tensile-loaded coupons , More work is required to determine if this is 
true for other configurations and loadings. However, if strain energy 
release rates must be separated into Gj, components using 

a numerical analysis, then the dependence of these calculated components 
on grid size should be carefiilly checked and documented. 

Immediate Results 

A simple rule of mixtirres analysis, using laminated plate theory, 
indicated that laminate stiffness was a linear function of delamination 
size. The analysis acciirately predicted stiffness loss due to edge 
delaminations in [±30/±30/90/90]„ graphite epoxy laminates. The linear 
stiffness relationship was used to derive a closed-form equation for the 
strain energy release rate, G, associated with delamination growth in 
unnotched laminates. The simple G equation was used to predict delamina- 
tion onset in ~ 1»2,3) laminates using a 

critical G determined from [±30/±30/90/90]o laminates. Stable 
c ° 

delamination growth in the [±30/±30/90/90]g laminates was characterized 



by developing a delamination resistance curve (R-curve). Delamination 
growth in [±30/±30/90/90]g laminates in fatigue was chajracterized by 
developing a power law correlation between G and delamination growth rates. 

Potential Applications 

Preliminary predictions of delamination onset in [+45v,/-^5„/0 /90„] 

laminates using a critical G„ determined from tests on [±30/±30/90/90]„ 

laminates indicated that G^ may be independent of the ply orientations 

that make up the delaminating interface. If this is true for stable 

delamination growth with increased tensile load and delamination growth 

in fatigue, then the delamination resistance curve (R-curve) and power 

law developed on [±30/±30/90/90] laminates can be used to predict 

s 

delamination growth in other laminates. 
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FINITE ELEMENT ANALYSIS 

Formulation 

The q.uasi -three-dimensional finite element analysis was developed 
in reference [7]. A displacement field of the form 

u = 0 qX + U(y,z) 

V = V(y,z) (Ai) 

w = W(y,z) 


was assumed, where was a prescribed uniform axial strain. Eight-noded 

quadrilateral, isoparametric elements with three degrees of freedom per node 
were used to model a cross section along the specimen length (fig. 18). Only 
one quarter of the cross section was modeled due to symmetry conditions. Each 
ply was modeled with one element through its thickness except for the cen- 
tral 90 ply, which was modeled with one element through its half-thickness. 
The graphite epoxy unidirectional properties [ll] used in the analysis were 

\l ~ ^38 GPa (20.0 Msi) 

®22 ^33 = 15 GPa (2.1 Msi) 

(A2) 

*^12 “ *^13 “ *^23 “ 5.9 GPa (O.85 Msi) 

^^12 ~ ^13 ~ ^23 ~ 0.21 
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Virtual Crack Extension Technique 

Previously, a virtual crack extension techniq.ue has been applied to 
finite element analysis of delaminations [l,4]. In this technique, the work 
required to close the delamination, expressed in terms of nodal forces and 
displacements, is assumed to be equivalent to the strain energy released as 
the delamination extends due to a constant nominal strain. This technique 
greatly simplifies the computation of Gj, Gjj, and Gjjj because 

knowledge of the singular stress field near the crack tip is not required, 

Figure 19 illustrates the technique as it was used with the quasi-three- 
dimensional finite element analysis. First, the nodal forces were calculated 
for an initial delamination of size a. Then, the delamination was extended 
an amount Aa and the resulting nodal displacements at the same location 
were calculated. The expressions for Gj, Gjj, and Gjjj in terms of 

npdal forces and displacements were 





■ \) * - «e) 

- ''c) " - v) 

- V) ♦ - 4 


(A3) 


where, for example, represents the force in the z direction at pode b 

b 

calculated from element 1. 

Figure 11 shows Gj and as functions of delamination size, 

G^^^ was negligible for this case. Strain energy release rates, G 

FEM 
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calculated using equations (A3), where normalized by the constant value pre- 
dicted from equation (lO). The total G, represented by Gj plus 
reached the value* predicted by equation (lO) once the deiamination had grown 
a very small distance in from the edge. Furthermore, the total G calculated 
was not sensitive to mesh refinement. However, the values of Gj and Gjj 
calciolated were sensitive to mesh refinement. 

Figure 11 shows the nodal discretization used in the finite element 
analysis. Four different initial delamination sizes, a, were modeled. For 
each a, the delamination was grown in ten increments, Aa, equal to one tenth 


*Because the finite element analysis models only one quarter of the 
laminate cross section due to symmetry, the analysis assumes four delaminations 
grow simultaneously, two on each side, located in the -30/90 interfaces. How- 
ever, in the derivation of equation (lO), only a single delamination was con- 
sidered on either edge to be consistant with the physically observed behavior 
shown in figure 3(c) and illustrated in figure 20(a). Making the same assump- 
tions as the finite element analysis would have yielded 


1/ = 2bilt 
A* = 4bJl 
A = hZa. 
dA = 4JI da 

which would result in G = ^^LAM " ®*)‘ Hence, the values calculated from 
the finite element analysis using equations (A3) were doubled to be consistent 


with equation (lO) and the observed physical behavior. 
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of the initial delamination size. As a was increased) the mesh refinement 
changed accordingly. For the smallest delamination size, continually 

increased as the delamination grew, but G^ was constant for all ten incre- 
ments of growth. The ten values of Gj for the smallest a are represented 
by a horizontal line segment. For the next three delamination sizes, both G^ 
and Gjj were constant for a^.! ten increments of growth, but the 
ratios Gj/Gjj changed. For the .largest delamination size modeled, Gj 
und Gjj were identical, both equal to one-half the total G. Hence, 
calculations of Gj and Gjj varied with mesh size. 
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EFFECT OF COUPLING ON STIFFNESS LOSS 
The stiffness of an arbitrary composite leuninate may be calculated from 
laminate theory [9 >10] as 


F - 1 

^LAM - XTHt 
11 

where 


(Bl) 


[A'] = [A]-^ + [A]-1[B][D*]-1[B][A]-1 
[D*] = [D] - [B][A]-1[B] . 

and [A], [B], [D] are the extensional, coupling, and bending stiffness 

matrices, respectively, defined in reference [9]. If the laminate is 
symmetric, the ,[b] matrix vanishes, and equation (B1) reduces to 
equation (l). 

For an arbitrary laminate containing delaminations in one or more 
interfaces, equation (Bl) may be used in equation ( 2 ) to calculate the 
stiffnesses, Ej^, of the sublaminates that are formed. For the 
[± 30 / 130 / 90 / 90 ]g laminate, with delaminations modeled in both - 30/90 inter- 
faces, equation (2) becomes equation (3) 

®^(±30) 2E(pq) 

E* = ^ i 

11 
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If th.e (^30)2 suDlaminate stiffnesses are calculated using equation (bi) 
and material properties frpm reference [ll], then E* = O.69 which 

exceeds the value of E* extrapolated from a least-squares regression line 
for the stiffness data. However, both -30/90 interfaces did not delaminate 
complete (fig. 3). As shown in figure 20 (a), delaminations shifted fi^om 
one -30/90 interface to another through various 90° ply cracks. Delaminations 
were believed to grow in this manner to reduce the effect of bending-extension 
cpupling that would have been present had only one interface been cleanly 
delaminateci or had both interfaces delaminated at the same time. Therefore, 
delaminations were modeled over both -3O/9O interfaces (fig. 20 (b)) but the 
binding-extension coupling in the (±30)2 sublarainates was negfected. 
fience , the stiffness of each sublaminate was evaluated using equation (l). 

This resulted in E* = 0.7^3 Ej^ which agreed with the value of E* 
extrapolated from a least-squres regression line for the stiffness data 
(fig. 10). 

In addition, stiffness loss for delaminations of various size was 
calculated from the finite element analysis. Concurrent strip delaminations 
were assumed in both -30/90 interfaces. Uniform axial extension was also 
assumed and, hence, bending-extension coupling was ignored. The stiffness 
corresponding to a particular delamination size was calculated using the 
technique outlined in reference [16]. Figure 21 shows the stiffness loss 
predicted by the linear rule of mixutes (R.O.M.) equation without coupling 
(E^ calculated from eq. (l)) and with coupling (E^ calculated from eq. (Bl)). 
The finite element predictions agreed with the rule of mixtures equation tb^at 
peglected bending-extension coupling. Hence, the finite element results also 
agreed with the da,ta plotted in figure 10. 
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STIFFNESS CHANGE FROM SECONDARY MECHANISMS 
Ninety Degree Ply Cracks 

The spacings of 90° ply cracks in the laminated portion of partially 
delaminated [±30/±30/90/90]g laminates were measured directly from enlarged 
x-ray photographs. The mean spacing is plotted in figure 22 as a function 
of delamination size (see fig. 9) for three quasi-static test specimens. As 
the delaminations grew, the crack spacings decreased to within values pre- 
dicted by a simple one-dimensional model [6]. 

Near the 90° ply cracks, load was assumed to be transferred to the 
neighboring plies and then back into the 90° plies to form the next crack. 

The model predicted that 90° ply cracks would reach a characteristic satura- 
tion spacing, K, determined by the laminate stacking sequence, ply thickness, 
and ply stiffness, but independent of load history. 

For the [±30/±30/90/90]g laminate, the model predicted a saturation 
spacing of 0.838 mm (0.033 in.). For a [±30/±30/902] sublaminate, repre- 
sentative of the delaminated portion of the [±30/±30/90/90]_ laminate, the 
model predicted a saturation spacing of 0,955 mm (O.O367 in.). 

Stiffness change due to 90*^ ply cracking in the partially delaminated 
specimen was estimated by reducing the maximum axial stress in the 90^ plies 
based on the density of cracks present. As shown in figure 23, the transverse 
stiffness of the 90^ plies, ^22 5 'was reduced by the percentage difference 
in the maximum axial stress, calculated from the one dimensional model, at 
the predicted crack spacing K and at one half the spacing K/2. This 
difference represents the reduction in load carried by the 90^ ply once the 
saturation crack pattern has formed. 



APPENDIX C 


For the [±30/130/90/90]^ laminate, a twelve percent reduction in the 
axial stress in the 90 plies was predicted. Reducing, the transverse 
stiffness, E 22 , of the 90 ° plies by twelve , percent resulted in a 0.7 percent 
reduction in [±30/±30/90/90]g laminate stiffness. Similarly for the 
[+ 30 / 130 / 902 ] sublaminate, a l4 percent reduction in the axial stress in 
the 9 o 9 plies was, predicted. Reducing the transverse stiffness, Egg, 
of the 90 ° pliep by l4 percent resulted in a 2.2 percent reduction in 
the [I 3 O/ 13 O/ 9 O 2 ] sublaminate stiffness. 

Because the estimated contributions of 90° ply cracks to stiffness 
reduction in the [l30/130/90/90]g delamination specimen were so small, 
they were neglected in calculating E*. However, for other laminates and 
materials, ply cracking may have a greater effect on axial stiffness and 
other stiffness parameters [lT,lS], 

+ 3 O/- 3 O delaminations 

The largest + 3 O/- 3 O interface delaminations, a, recorded during the 
four quasi-static tension tests were measured and tabulated in table 2 . 

Only specijnen E20 had a + 30/-30 delamination of significant size. Also 
shown in table 2 are predicted stiffness values for the [l 30 /± 30 / 90 /^] 

s 

laminate j containi^ng both -30/90 and +30/— 30 interface delaminations. The 
stiffness was predicted using equations (2) and (4) first for the 
+ 3 O/- 3 O delamination, and then for the -30/90 delamination. The stiffness 
predicted for specimen E20 was three percent less than the stiffness pre- 
dicted assuming only -30/90 delaraination. However, the data agreed best 
with the prediction that omitted the +30/-30 delamination* Hence, the 
effect of + 3 O/- 3 O interface delaminations on measured [±30/±30/90/9^] 

s 

stiffness loss was neglected in calculating E*. 
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TABLE 1.- FATIGUE GROWTH RATE MEASUREMENTS 


^max 
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■'max 

. 2 ^ 


Spec . 

# 


da' 


mm 


dN* cycle 


Area 

measvire 


Stiffness 
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dN’ 


cycle 


Cycle range of 
linear fit 
(cycles X 10^) 


30 ( 

DO 

1 C 
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22.3 
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F2h 
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\ 
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80 
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51.2 

- 40.3 
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r 
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46.9 




3500 , 

138 
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52.6 
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87.1 
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64.6 
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47.9 

46.6 
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TABLE 2.- EPPECT OP +30/-30 DELAMINATION ON STIPPNESS PREDICTION 


Specimen 

# 

-30/90 

+30/-30 


prediction 

E/Eq measured 

a 

b 

a 

b 

■With S i 

t 

Without a 

P9 

0.81^8 

0.009 

0.778 

0.780 

0.782 

E20 

00 

o 

oo 

• 

.109 

.768 

.792 

1 .806 

C18 

.812 

.OOlt 

.790 

.790 

.798 

D6 

.773 

O 

o 

.800 

0 

0 

OO 

.787 
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Figure 2.- Radiographs showing damage growth through the width of [±30/±30/90/90 ] 5 

laminates. 
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Damage location through the t 






TERFACES 

MATED 


+30/-30 interfaces 



0 ^, GPa 


Figure 5.- Through-thickness distribution at edge of 

[±30/±30/90/90 ] 3 laminate. 
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Figure 6.- Through-widtti distributions of and a, near the edge. 
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Figure 7.- Rule of mixtures analysis of stiffness loss. 
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Figure 9.- Illustration of strip delamination approximation. 
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Figure 11.- Strain energy release rate components from finite element analysis as a 

function of delamination size. 
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Figure 13.- Delamination resistance curve (R-eurve) for [±30/±3O/9O/9O3g T300-5208 

graphite-epoxy laminates. 
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Figure 17." Effect of delaminotion measurement method on power law curve fit. 




Figure 18.- Finite element discretization of specimen cross section. 
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Figure 19.- Virtual crack extension technique used to calculate 

strain energy release rates. 
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Figure 21.- Effect of coupling on -stiffness loss prediction. 
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Figure 22,- Ninety degree ply crack spacing as n function of 

delaniination size. 
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Figure 23.- Estimate of stiffness loss from 90° ply cracks. 
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